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Within the QCD factorization formalism, we study the possible impacts of the nonuniversal Z' 
model, which provides a flavor-changing neutral current at the tree level, on rare decays B — > Kq n. 
Under two different scenarios (SI and S2) for identifying the scalar meson X"q(1430), the branching 
ratios, CP asymmetries, and isospin asymmetries are calculated in both the standard model (SM) and 
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the family nonuniversal Z' model. We find that the branching ratios and CP asymmetries are sensitive 



to weak annihilation. In the SM, with Pa = 1 and <j>A G [—30°, 30°], the branching ratios of SI (S2) are 
smaller (larger) than the experimental data. Adding the connibution of the Z' boson in two different 

m : 

cases (Case-I and Case-II), for SI, the branching ratios are still far away from experiment. For S2, 
' in Case-II, the branching ratios become smaller and can accommodate the data; in Case-I, although 

the center values are enhanced, they can also explain the data with large uncertainties. Similar 
conclusions are also reached for CP asymmetries. Our results indicate that S2 is more favored than 
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SI, even after considering new physics effects. Moreover, if there exists a nonuniversal Z' boson, 
Case-II is preferred. All results can be tested in the LHC-b experiment and forthcoming super-B 



I. INTRODUCTION 



Recently, with rich events in two B factories, measurements of B meson nonleptonic charmless decays 
involving scalar mesons have become available. Among these decays, the processes B — ^ K^n are attractive 
since they are dominantly induced by the flavor-changing neutral current (FCNC) transition b — > sqq {q = 
u,d,s). Such a transition forbidden at the tree level in the Standard Model (SM) is expected to be an 
excellent ground for testing SM and searching for new physics (NP) beyond SM. Therefore, many similar 
decay modes induced by FCNC have been explored widely in the literatures, such as B — > KK,Kr]^' ,$K(*\ 
The recent reviews can be found, for example, in Ref. [1]. For the concerned decay modes B — > K^n, the 
latest world averaged branching ratios from Heavy Flavor Average Group (2|] are listed as: 

BR(B + ^^°( 1430) 71+) = (45.1 ±6.3) x 10" 6 ; 
BR{B°^K* + { 1430)71-) = (33.5+^) x 10~ 6 ; 

B/?(B°^^°(1430)7r ) = (11.7±f;|) x HT 6 . (1) 
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Direct CP asymmetries of above decays have also been measured recently by BaBar and Belle experiments, 
which will be shown in Sec.|lVl As direct CP violation is sensitive to the strong phase involved in the decay 
process, the comparison between theory and experiment will offer us information on the strong phases 
necessary for producing the measured direct CP asymmetries. Comparing the predicted results of the SM yO 
with experimental data, ie. Eq.(Q~|), we notice that the theoretical results cannot accommodate the data well 
even with large uncertainties. So, it is worth while to explore whether some new physics models could 
explain the data. 

When discussing the B meson non-leptonic charmless decays, the hadronic matrix elements are required. 
In the past few years, several novel methods have been proposed to study matrix elements related to exclu- 
sive hadronic B decays, such as naive factorization (NF)J4|], generalized factorization (3], the perturbative 
QCD method (pQCD) |fj], QCD factorization (QCDF) ffl\, the soft collinear effective theory (SCET) |8D, 
and so on. Among these approaches, QCDF based on collinear factorization is a systematic framework to 
compute these matrix elements from QCD theory, and it holds in the heavy quark limit mi, — > °° and the 
heavy quark symmetry. Thus, we shall use QCDF approach in the following calculations. 

Although the study of scalar meson spectrum has been an interesting topic for a long time, the underlying 
structure of the light scalar meson is still controversial until now. In the literature, there are many schemes 
for the classification of them. Here we present two typical scenarios to describe the scalar mesons (9J]. 
Scenario- 1 (SI) is the naive 2-quark model: the nonet mesons below 1 GeV are treated as the lowest lying 
states, and the ones near 1.5 GeV are the first orbitally excited states. In scenario-2 (S2), the nonet mesons 
near 1.5 GeV are regarded as the lowest lying states, while the mesons below 1 GeV may be viewed as 
exotic states beyond the two-quark model. Since the mass of ^,(1430) is very near 1.5 GeV, thus it should 
be composed by two quarks in both SI and S2, but the decay constants and distribution amplitudes are 
different in the different scenarios. Under abovepictures, the two body nonleptonic B decays involving 
scalar mesons have been explored in both QCDF olljll H] and pQCD approaches 1 12- l3l. 

As stated before, B — > KqK decays are dominantly induced by FCNC b — >• sqq transition, hence they 
are sensitive to new physics contributions even if they are suppressed by a large mass parameter which 
characterizes the new physics scale. To search for signals of NP, a model independent analysis is not 
suitable for the current status. It is the purpose of this work to show that a new physics effect of similar size 
can be obtained from some models with an extra Z' boson. Z' bosons are known to naturally exist in certain 



well-motivated extensions of the SM, such as the string theory [18], the grand unified theories [ 19], the little 
Higgs models[20], light U-boson model [2l| . by adding additional U{\)' gauge symmetry. Among those 
models, a well-motivated Z' model for low energy systems is the so-called family non-universal Z' model, 
where the Z' couplings are affected by fermion mixing and are not diagonal in the mass basis. Non- trivial 
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flavor-changing neutral current (FCNC) effects at the tree level mediated by the Z' therefore are induced, 
which play an important role in explaining the CP asymmetries in the current high energy experiments 
by introdu cing new weak phases. The effects of Z' boson in B sector have been investigated in a number 



of papers [22 



26(1 . In this work, we will show the implications of the family nonuniversal Z' model on 



B — > KqTI decays. 

The layout of this paper is as follows. In Secfll] we firstly present the formulaes of B — > Kq % in the 
SM within the QCDF approach, involving the effective Hamiltonian and the amplitudes. In SecfnTJ we 
specify our flavor-changing Z' model, and how the effective Hamiltonian responsible for hadronic B decays 
is modified. The numerical results and discussions are given in SectJIVIThe conclusions are presented in 
the final section. 



II. CALCULATION IN THE STANDARD MODEL 

In the two-quark picture of SI and S2, the two kinds of decay constants of scalar meson S are defined 

by: 

(%)|? 2 y M ?l|0) = fsPfi, (S(p)\q 2 qi\0)=m s f s . (2) 

The vector decay constant f$ and the scale-dependent scalar decay constant f$ are related by equations of 
motion 

rat 

Ms/s = fs, with ju 5 = — 7-v 7-t, (3) 

m2(H)-mi(n) 

where 1112 and m\ are the running current quark masses. Therefore, contrary to the case of pseudoscalar one, 
the vector decay constant of the scalar meson, namely, fs, will vanish in the SU(3) limit. In other words, 
the vector decay constant of ^(1430) is fairly small. 

As for the scalar meson wave function, the twist-2 and twist-3 light-cone distribution amplitudes 
(LCDAs) for different components could be combined into a single matrix element: 

<**+ (p) | up (z)s a (0) |0> = -J= jf dxJ** (x) + m s ^ + (x) + -msG^p^^ (x) | . (4) 



The distribution amplitudes (j>iq(x), an d ^k^ x ) are normalized as: 







1 dxfa* (x) = I' dx& (x) = [ dxQZ, (x) = ^=, (5) 



2^6' Jo K ° Jo K ° 2^6' 
and (jc) = The twist " 2 LCDA 

can be expanded in the Gegenbauer polynomials: 

1 CO 

te(x,n) = -=f s (n)6x(l -x) £ B m {ji)&d\2x-\). (6) 



c m=\ 
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The decay constants and the Gegenbauer moments for twist-2 wave function in two different scenarios 
have been studied explicitly in Ref. I3I [k]] using the QCD sum rule approach. As for the explicit form 
of the Gegenbauer moments for the twist-3 wave functions, there exist few drawbacks in the theoretical 
calculation J27II . thus we choice the asymptotic form for simplicity: 



0j 



1 



=//, 4>s 



1 



=/s(l-2x). 



For the pion meson, the asymptotic forms for twist-2 and twist-3 distribution amplitudes are also adopted 



(7) 



0p (*) = fpte( 1-x), (j)p (x) =f P , (f>p (x) = f P 6x( l-x). 



(8) 



The form factors of B — > P, S transitions are defined by [4J] : 



(p(p'm\B( P ))= p. 



m% — nit 
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7 BSfJ2\ 1 m B~ m S„ pBS{2\ 



(9) 



T J T 

where P^ = (p+p')^, q^ = (p— p')u . Various form factors have been evaluated by utilizing the relativistic 
covariant light-front quark model 12811 . And the momentum dependence is fitted to a 3-parameter form 

F(0) 



F{q Z ) 



(10) 



The parameters a and b relevant for our purposes are refereed to Ref. j28ll . 

Although we concentrate on the study of new physics, the used notation for new interacting operators 
will be similar to those presented in the SM. Therefore, it is useful to introduce the effective operators of 



the SM. Thus, we describe the effective Hamiltonian for b — > sqq decays as 



V p=u,c 



C 1 ( At )0 ( / ) ( At )+C 2 ( At )Of(At) + lQ(M)0 ( -(M) 



(<?), 
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1=3 



(11) 



where X q = VghV L are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and the operators 0\-0\q 
are defined as 1 29D 

( i ] = {saqp)v-A(q~i}ba)v-A , o!f = {s a q a )v-A{qpbp) V - A , 
O3 = {s a b a )v-AY^(^li1p)v-A , O4 = (s a bp) V - A Y,(qpq<x)v-A , 

1 9 

O5 = (s a b a )v-AY,(^p9p)v+A , 6 = (sabpjv-A^iMpqcdv+A , 



1 1 

3 3 

Oi = ^(s a b a )v-AY, e q(^plp)v+A , 08 = y(s a bp) V -AY, e q(^pla)v+A , 

9 1 

3 3 

O9 = ^{s a b a )v~AY J e 'Mp c lp)v-A , 010 = ~(Sabp)v-AY, e <l(QpQct)v-A , 
1 q 1 q 



(12) 
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with a and /3 being the color indices. In Eq.dTTb. 0\-02 are from the tree level of weak interactions, 03-06 
are the so-called QCD penguin operators and Oj-Oio are the electroweak penguin operators, while C\-C\q 
are the corresponding Wilson coefficients. 

In the QCDF approach, the contribution of the non-perturbative sector is dominated by the form factors 
and the non-factorizable impact in the hadronic matrix elements is controlled by hard gluon exchange. The 
hadronic matrix elements of the decay can be written as 

{M x M 2 \Oi\B) = £ F ?^ Ml jf dxTj(x)* Mi (x) + jf d£, £ dx jf dyT/ 1 (t;,x,y)<S> B (t; )$ M] (x)<D M2 (y). (13) 

Here Tf, and T? 1 denote the perturbative short-distance interactions and can be calculated perturbatively. 
<J>x (x) are non-perturbative light-cone distribution amplitudes, which should be universal. Using the weak 
effective Hamiltonian given by Eq.(fTT1) and the definitions of a, and bi in Ref . O, lZ[| , we can now write the 
decay amplitudes of B — > KqK as: 

1 



A{B~ 



V z p=u,c IV z / nic 



f K *F B7l (m 2 K M B -mi) 



+fB(b 2 S^ + b 3 +b 3fiW ) >, 



(14) 



A(B~^K^K U ) 



Gf 
2 



m B -m n ) 



AF ( f ^ (m 2 J (ml - ,1^ ) + f B (b 2 8£ + b 3 + Z? 3 ,ew 



a 2 8? + -{a 9 -a 7 



(15) 



V 2 p=MjC [ v 7 



/fi (*3 - ~^3,EW 



(16) 



z p=«,c IV z / jtx; 



«2<5,f + -(09-07) 



fn F o ( m l)( m l ~ mj~ )+f B {-h + ^h,Ew] 



71KZ 



(17) 



where X p = V p bV* and 



2rai» 



m b (n)(m s (n)-m q (n))' 



(18) 
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In the above formulaes, the order of the arguments of the af (M1M2) and bi{M\M-i) coefficients is dictated 
by the subscript M1M2, where is the emitted meson and M\ shares the same spectator quark with the 
B meson. For the annihilation diagram, Mi is referred to the one containing an anti-quark from the weak 
vertex, while M2 contains a quark from the weak vertex. Note that the coefficients a,- come from vertex 
corrections and hard spectator corrections, and bj represent of contribution of annihilation diagrams. Both 
a, and bj can be found in Ref.[3]. It must be emphasized that we shall evaluate the vertex corrections to 
the decay amplitudes at the scale /J. = nit,/ 2. In contrast, the hard spectator and annihilation contributions 
should be evaluated at the hard-collinear scale /I/, = ^/J^Kh with A/ t « 500 MeV. 

In QCDF approach, the annihilation amplitude has endpoint divergences even at twist-2 level and the 
hard spectator scattering diagram at twist-3 order is power suppressed and posses soft and collinear di- 
vergences arising from the soft spectator quark. Since the treatment of endpoint divergences is model 
dependent, subleading power corrections generally can be studied only in a phenomenological way. We 
shall follow lO, 3] to parameterize the endpoint divergence = Jq dxjx in the annihilation diagram as 

* A=in G?) (i+p ^ A) ' (19) 

with the unknown real parameters and fa. Likewise, the endpoint divergence Xh in the hard spectator 
contributions can be parameterized in a similar manner. In the SeclTVl we will see that such divergence is 
the main source of the uncertainty for the concerned decay modes. 



III. THE FAMILY NON-UNIVERSAL Z' MODEL 



As mentioned before, a family non-universal Z' model leads to FCNC at the tree level due to the non- 
diagonal chiral coupling matrix, which makes itself become interesting in some penguin dominate pro- 
cesses. The basic formalism of flavor changing effects in the Z' model with family nonuniversal and/or 



nondiagonal couplings has been laid out in Refs. [2: 



3 



2611 . to which we refer readers for detail. The detailed 



phenomenological analysis for various low energy physics, especially for B meson decays, could be found 



in Refs. (23 



- 25I1 . Here we just briefly review the ingredients needed in this paper. 
In practice, neglecting the renormalization group (RG) running between m w and m z > and mixing between 
Z' and Z boson of the SM, we write the Z' term of the neutral-current Lagrangian in the gauge basis as 



J? = -g>J>Z>», 



(20) 



where g' is the gauge coupling constant of extra U(l)' group at the electro- weak mw scale. The chiral 
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current J' is expressed as: 



(Bf,) J P L + (B«.)^ 



Yj, (21) 



where the chirality projection operators are Plr = (1 ± Ys)/2 and Bfj refers to the effective Z' couplings 
to the quarks i and j at the electroweak scale. For simplicity, we assume that the right hand couplings are 
flavor-diagonal and neglect Bf b . Compared with Eq.dTTI). the effective Hamiltonian for b — > sqq transition 
with Z' boson can be written as 

^eI' = ^(^) 2 ^(^%-Al[^(w)y-A + ^(w)y+A]+h.c., (22) 

where mz> is the mass of the new gauge boson. In fact, the forms of four-quark operators in Eq. (1221 already 
exist in the SM, so we rewrite it as 

J0>£ = -^Lv tb V*Y(AC 3 O q 3 +AC 5 Ol +AC 7 O q 7 +AC 9 Ol) +h.c, (23) 
V2 q 

where Of (i = 3,5,7,9) are the effective four-quark operators in the SM. AC, denote the modifications to 
the corresponding SM Wilson coefficients, which are expressed as 



W tb V t * y gl m z , 



A C 97 = -JBtf-Bff), (24) 

Generally, the diagonal elements of the effective coupling matrices B qq are expected to be real as a con- 
sequence of the hermiticity of the effective weak Hamiltonian. However, the off-diagonal one B sb perhaps 
contains a new weak phase We also suppose B L qq = B qq = B qq , so as to reduce the new parameters. For 
convenience we can represent AC/ as 

AC 3 , 5 = 2^£e*, AC 9J = 4^f^*, (25) 
VtbVts V t bV ts 

where £ and t, are defined, respectively, as 

It is stressed that the other SM Wilson coefficients may also receive contributions from the Z' boson through 
renormalization group (RG) evolution. With our assumption that no significant RG running effect between 
M' z and M\y scales, the RG evolution of the modified Wilson coefficients is exactly the same as the ones in 
the SM 1291]. 
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In order to show the effects of Z' boson clearly, our analysis are divided into the two cases with two 
different simplifications, 

' B m = -2B dd , C = 0, § =X, Case I; 
B m =B dd , £ = -X, §=0, Case II. 



(27) 



with 



B sh B dd 



VtbVu 



■y 



B L sb B dd 



VtbV t * 



(28) 



Thus, there are only two parameters, X and weak phase (j) s left, in the sequential numerical calculations and 
discussions. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

To obtain the numerical results, we list the parameters related to the SM firstly. As stated in Section. H 
because we have not a clear conclusion whether Kq (1430) belongs to the first orbitally excited state (SI) or 
the low lying state (S2), we have to calculate the processes under both scenarios. So, the decay constants, 
Gegenbauer moments, and form factors in different scenarios are listed as follows (3D : 

SI: A; (1.0 GeV) = -300 MeV;/^(2.1 GeV) = -370 MeV;fli(1.0 GeV) = 0.58; 
fli(2.1 GeV) = 0.39; B3 (1.0 GeV) = -1.20; B 3 (2.1 GeV) = -0.70; F** 3 (0) = Ff *"(()) =0.21; (29) 

S2: fg* ( i430) (1.0 GeV) =445 MeV;/^ (1 4 30) (2.1 GeV) = 550 MeV;*i(1.0 GeV) = -0.57; 
B x (2.\ GeV) = -0.39;fl 3 (l-0 GeV) = -0.42; B 3 (2.1 GeV) = -0.25; F<f*°(0) = if*°(0) = 0.26(30) 

Now that the uncertainties for the above parameters have been explored explicitly in Ref.[3], and we will 
not discuss the errors caused by them in the current work. 

In Ref. [3], the authors concluded that the theoretical errors are dominated by the l/mj power corrections 
due to the weak annihilations. Moreover, the weak annihilation contributions to B — > SP could be much 
larger than the B — > PP case, because the helicity suppression appeared in the B — > PP case can be alleviated 
in the scalar production with the non-vanishing orbital angular momentum in the scalar state. In order to 
accommodate the data, one has to take into account the power corrections due to the p# and p A from the 
hard spectator interactions and weak annihilations, respectively. In Ref. [3], Cheng et.al found that the 
predictions are far away from the experimental data if by setting = 0, which indicates that will be 
nonzero. Meanwhile, for B — > PP,PV modes [7], the errors due to weak annihilations are comparable to or 
much smaller than the center values, and the fitting results show that Pa = 1 and A = 0°. Hence, in this 
work, we adopt Ph = Pa = 1, and set the strong phases in the ranges [—30°, 30°]. 
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TABLE I: Branching ratios (in units of 10 ) of B — > Kq k in the SM and the non-universal Z' model. 



Decay Mode 



SM 



SI 

Case I 



Case II 



SM 



S2 
Case I 



Case II 



Expt 



B 
B 

fi° - 



— *o 



r v 



23.0: 


H.2 
-5.9 


25.7 


+5.0+2.8 
-4.7-7.7 


17.21 


9.31 


1.0 
1.9 


17.9 


1-0.8+11.4 
-3.3-17.3 


6.81 


21.3; 


1-0.7 
-5.1 


27.2; 


f0.2+6.8 
-6.4-15.6 


15.91 


12.9; 


(-0.3 
-3.6 


9.4j 


0.3+12.0 
2.7-2.9 


9.81 



+ 1.3+19.6 



4.4 
0.3-1 



-5.0 

-8.7 



-18.4 
-4.6 

10.3 



74.7; 
38.9 
70.0 
33.6 



-1.0 
20.6 

f0.4 
-8.9 

-0.6 
-17.2 

0.4 

-9.8 



93.8 
75.3 



-1.7 

-25.: 



-20.6 
-50.6 



-0.0+46.8 
-14.5-73.1 



S o O+0.0+13.4 
OJ - : '-16.2-36.3 



53.81 

28.2 



50.2 



1.3+16.5 

17.5- 17.4 

+0.5+8.4 
-7.7-8.9 

+0.9+15.6 

14.6- 16.4 



45.1 



h6.3 
-6.3 



22.0 



-0.5- 
-7.2- 



-41.9 
8.5 



23.8 



-0.5+7.7 
-8.2-8.0 



33.5 
11.7 



-3.9 
-3.8 



h4.2 
-3.8 



TABLE II: CP asymmetry (in %) of B — > KZ % in the SM and the non-universal Z' model. 



Decay Mode 



SM 



SI 
Case I 



Case II 



SM 



S2 
Case I 



Case II 



Expt 



5°- 



■K*-7t° 



K*-n+ 

r v 



1.0 

-0.5 
2.0 
3.1 



1.9 
1.9 

+3.8 
-2.6 

+2.7 
-3.9 

-2.6 
-4.0 



1.0 



+ 1.4+0.2 
1.6-0.1 



h2.4- 
-2.1- 



0.2 
0.4 



-0.4 



1.7 

3.7 



-2.7+0.6 
-2.0-2.8 

+2.3 + 1.5 
-2.6-0.7 



1.2; 

— 6 1 : ,M 
2.4 



2.8-0.2 
+3.1 + 1.0 



f -4.6-1.4 



0.06+g;« 



-2.4 
-2.9 

+2.9 
-2.7 



06 +0 - 4+0 - 3 

U - UO -0.5-0.3 



0.03 



h0.66- 
-0.82- 



0.03 
0.07 



1.0 
-0.8 



H3.0H 
-4.7- 



1.3 

2.0 



3.7 



h2.9H 
-4.5- 



1.0 
1.8 



i q+4.0 



0.7 



-0.7 
-2.5 



+ 1.8+5.4 
-2.0-1.2 

+2.1+0.6 
-2.2-0.9 

-5.3+3.3 
-4.4-2.2 



,+2.7+0.3 
--3.5-0.2 



1.1 



+3.6+0.2 
-3.3-0.5 



-2.5 H 



0.5 
1.0 



-7 
-34 



14 

-14 

19 
19 



With above parameters, we present our predictions of the SM in Tablefl] under two different scenarios. 
For the center values, we also assign (j) A = (j) H = 0. In order to obtain the errors, we scan randomly the points 
in the ranges § A G [-30°, 30°] and (j> H G [-30°, 30°]. So, the only theoretical errors of the SM results are 
due to the strong phases 0^ and <pH. Because we fully consider the weak annihilations, our results are much 
larger than those in Ref.[3], especially for the center values. Compared with the data, the theoretical results 
in this work are still much smaller (larger) than the data under two scenarios, except for mode 6° — > K^n . 
If one want to fit the data absolutely, p A ~ 1.3 for SI and p A ~ 0.7 for S2 are required, respectively, which 
are a bit larger/smaller by 30% than the fitted results from B — > PP,PV. Compared with predictions of 
Ref. [15] obtained in the pQCD approach based on kj factorization, our results are a bit larger than theirs 
in S2, but agree with their results in SI with large uncertainties. 

We next turn to the implications of the non-universal Z' model for the B ^ K^n decays. Let us firstly 
consider the range of X, which is the most important parameter in this model. Generally, we always expect 
g' /gi ~ 1, if both the U(l) gauge groups have the same origin from some grand unified theories. Mz/M z / ~ 
0. 1 for TeV scale neutral Z' boson is also expected so as to the Z' could be detected in the running Large 



Hadron Collider (LHC), which results in y ~ 10 2 . In the first paper of Ref. 12311 assuming a small mixing 
between Z — Z' bosons, the value of y is taken as y ~ 10~ 3 . In order to explain the mass difference of 
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B s —B s mixing, \B^ b \ ~ |^foV ; *| is required. Similarly, the CP asymmetries in B — > <pK, %K can be resolved 
if \B^ h Bss R \ ~ |V;/,V,*|, which indicates \B^ \ ~ 1. Above issues have been discussed widely in Ref. 12411 . 
Summing up above analysis, we thereby assume that X G (10~ 3 , 10~ 2 ). For weak phase ij> s , though many 
attempts have been done to constrain it 12511 . we here left it as a free parameter. 

The calculated results for branching ratios with two different cases in the family non-universal Z' model 
are also exhibited in TableJH and for the center values we use X = 0.005 and <j> s = 0°. To obtain the second 
errors, we also scan randomly the points in the ranges X £ [0.001,0.01] and (j> s E [—180°, 180°], while the 
first errors come from the weak annihilations. The table shows to us that the two cases of Z' models can 
change the branching ratios remarkably in the two different scenarios. It is clear that the Z' will enhance 
the branching ratios in Case-I, while in Case-II the branching ratios are decreased. The reason is that the 
variation tendencies of Wilson coefficients are different in the two different cases, which could be seen in 
Eq.(25) and (27) easily. 

For SI, the branching ratios of the first three decay modes cannot agree with data unless the upper limits 
in Case-II of the Z' model are taken. Unfortunately, with the upper limit values, the branching ratio of 
B — y Kq ti is much larger than the experimental data. For S2, the branching ratios with a Z' boson can 
accommodate experimental data well in two cases with large uncertainties. If we care about the center 
values very much, it seems that results of Case-II are preferable. If further theories and/or experiments can 
confirm the existence of Z', one could correspondingly cross-check the couplings and the mass of it with all 
above results in turn. 

In the experimental side, another important observable in B physics is CP asymmetry, in particular of 
the direct CP asymmetry. In TableuTl we list the direct CP asymmetries of concerned modes in different 
scenarios and different cases of the Z' model. Generally, the strong phases calculable in the QCD factor- 
ization are so small that the CP asymmetries are at most a few percent, as shown in the table. In SI we 
note that the center values have different signs with the experimental data. Adding the Z' contribution, 
although the large uncertainties perhaps alleviate the disparity of B — > K K~ , but the large asymmetries 
in the fi° K*-tc+ and fl° K*°7l° cannot be explained yet. In S2, for fl° K*-7C+ and fi° K*°7C°, 
the signs of center values are same as those of data. Furthermore, the CP asymmetries of B — > K n~ in 
the SM and Z' model are almost null, which are close to the upper limits of experiment. Considering the 
large uncertainties, the results of S2 in both SM and Z' models can accommodate the data, except for the 
unexpectedly large asymmetry of B° — > Kq°k°, which should be measured critically in future. However, as 
pointed out in Ref. [30], final state interaction may have important effects on the decay rates and their direct 
CP violations, especially for the latter. However, this is beyond the scope of the present work. 

Let us now analyze the impact of Z' on the isospin symmetry breaking. To explore the deviation from 
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TABLE III: Ratios of the branching fractions in the SM and the non-universal Z' model. 



SM 



SI 

Case I 



Case II 



S2 

SM Case I 



Case II 



Expt. 



Ri 

R 2 
R3 



0.60+°;°° 0.351 
0.40+°$ 0.70: 



0.03- 

0.02- 

0.04- 

0.09- 



1 nn+ 03 O 87 +0-28- 
i ' uu -0.02 u '°' -0.05-0.10 



1.50 
0.15 

0.33 
0.66 

1-0.59 



62 +0 °- 

u ' oz -0.05- 

n oq+0.05- 
u - J ^-0.00- 



0.02 
0.03 

0.02 
0.01 



1.00 



+0.04+0.01 
-0.02-0.01 



0.48 
0.52 
0.99 



hO.01 
-0.03 

-0.03 
-0.01 

-0.02 
-0.03 



n o/r+0.02+1.09 
U ' ZO -0.04-0.12 



-0.09+0.26 
-0.05-0.75 



rv A o+0.01+0.00 
U ' HO -0.03-0.00 



0.35 



-0.18 
-0.15 



0.80 
1.05 



+0.08+0.27 
-0.10-0.21 



0.52 
0.99 



+0.04- 
0.01- 

+0.02H 
-0.04-0.00 



0.01 
0.00 

0.01 



1.25 



-0.36 
-0.29 



the isospin limit, it is convenient to define the following three parameters: 

Zrc(5°->^ (1430)ft ) 

^1 = =n ' ^1) 

BR(W -»-Xj-(1430)w + ) 

BR(B-^K*-(l430)7t°) 

Ri = — — , (32) 

BR(B- ^Tq( 1430)71-) 

R3 = BR ( B ->^°(1430)k-) 

3 BR^^iq- (1430)jt+)' 

Because they are the ratios of the branching fractions, they should be less sensitive to the non-perturbative 
inputs than other observables discussed before, therefore it is more persuasive to test them in both theoretical 
and experimental sides. In the isospin limits, i.e., ignoring the electroweak penguins, R\ , R2 and /?3 are equal 
to 0.5, 0.5, and 1.0, respectively. So, the deviations reflect the magnitudes of the electroweak penguins 
directly. The results of SM and the non-universal Z' model are listed in TablelDTl In the SM, it appears that 
the deviations from the isospin limit are not large in both scenarios, which shows that the QCD penguins 
are dominant. For Case-I of the Z' model, the new physics just revise the Wilson coefficients of electroweak 
penguin operators, which could break the isospin symmetry. So, the ratios will be changed remarkably in 
both scenarios, as shown in the table. In Figured] we also present the variations R12 as functions of the 
new weak phase (j) with different X = 0.001,0.005,0.01 in SI (up panels) and S2 (down panels), so as to 
show the effect of two parameters X and (p. From the figures, we see that the R12 change remarkably when 
X = 0.01 and 0.005. As X = 0.001, 7?+2 almost have same values as predictions of the SM. For Case- 
II, the Z' boson changes the Wilson coefficients of QCD penguins, so the isospin symmetries are almost 
unchanged, as shown in TablejETl To sum up, the measurements of the Ri will help us determine whether 
QCD or electroweak interactions will be changed and then test the corresponding new physics models. 

Finally, we will go back to the discussion of two scenarios. As aforementioned, ^(1430) is regarded as 
two-quark state in both SI and S2, but the only controversy is whether it belongs to ground state or the first 
excited state. Through calculation and comparison above, we favor the second scenario, which means that 
^*(1430) is the lowest lying qq state. Namely, the scalar mesons lower than lGeV are four-quark states. 
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FIG. 1: R\ and Ri as functions of weak phase <j> s with different X in different scenarios and cases. 



This conclusion is also consistent with those of Refs.[3, 12, 



V. SUMMARY 



Based on the QCD factorization approach, we have investigated in this work B ^ K^k decays in the 
SM and a family non-universal Z' model. Because the inner structure of ^,(1430) is not clear enough, we 
calculated the branching ratios under two different scenarios (SI and S2). After calculation, we found that 
the branching ratios are sensitive to the weak annihilations. In the SM, with = 1 and <p A G [—30°, 30°], 
the branching ratios of SI (S2) are smaller (larger) than the experimental data. Considering the Z' boson in 
two different cases, for SI, the branching ratios are still far away from experiment. For S2, the branching 
ratios become smaller and can accommodate the data in Case-II; in Case-I, the results can also explain the 
data but with large uncertainties. Furthermore, the other interesting observables, such as CP asymmetries 
and isospin asymmetries, are also calculated. Compared with data, we favor that ^(1430) is the lowest 
lying qq state. Moreover, if there exists a Z' boson, Case-II is preferable. All above results will be tested in 
the B factories, LHC-b and the forthcoming super-B factory. 
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